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Abstract 

The air fluorescence technique has long been used to detect extensive air showers 
and to reconstruct its geometry and energy. The fluorescence photon yield of an 
electron in air is of main importance in the reconstruction procedure. Historically, 
the fluorescence yield used in the reconstruction of the showers is approximated at 
all energies by that of an electron with kinetic energy of 80 MeV, because these 
are the most abundant in a shower. In this paper, we calculate the fluorescence 
yield taking into account the energy spectrum of electrons in showers initiated by 
proton and iron nuclei as a function of height. We compare our results with previous 
calculations based on mono energetic electrons (80 MeV) and a difference in excess 
of 8% is found. Finally, the influence of a more realistic fluorescence yield in the 
shower energy reconstruction is also discussed. 
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1 Introduction 



Fluorescence telescopes have been successfully used to measure extensive air 
shower since the Fly's Eye [1] experiment proved its efficiency. Presently, the 
HiRes [2], Auger [3], EUSO [4], OWL [5] and Telescope Array [6]) Experiments 
are using or planning to use the fluorescence technique to study the cosmic 
ray with energy above 10^^ eV. 
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The detection of the shower by a fluorescence detector is performed by mea- 
suring the amount of fluorescence light produced in air by the excitation of 
the nitrogen molecules as particles (essentially electrons and positrons) travel 
along the atmosphere. The number of photons detected by the telescope can 
be converted to the number of charged particles in the shower at a given 
height taking into account the absorption and scattering of the fluorescence 
light along the path from the emission to the telescope and knowing the fluo- 
rescence yield of the particle. 

Prom the number of particles as a function of height one is able to reconstruct 
the energy of the shower following the arguments in reference [7]. Therefore, 
the fluorescence yield has a direct influence in the energy estimation of a 
shower. 

The importance of knowing the fluorescence yield of the electron as a func- 
tion of energy is underlined by the number of experiments presently under 
operation to measure the yield as can be seen in references [8,9,10,11]. 

Prom the measurements of the yield in the experiments cited above to the 
use of the fluorescence yield in cosmic rays showers a number of assumptions 
are made. Specially, the fluorescence yield of the particles in a shower are 
considered to be the yield of electrons with kinetic energy of 80 MeV. 

The contribution of other particles, like muons and hadrons and particles with 
energy below the simulation energy cut, to the total productions of fluores- 
cence light in a shower has been extensively studied in a recent paper by H. 
Barbosa et al. [12] . However, the mono energetic approximation has survived 
widely unchallenged. 

The scope of the present work is to evaluate the approximation of the energy 
spectrum of electron by single energy (80 MeV) in the evaluation of the fluo- 
rescence yield. We calculate the fluorescence yield taking into account the true 
energy spectrum of electrons (= electrons -|- positrons) in showers initiated by 
proton and iron nuclei with energy 10^^, 10^^ and 10^° eV. Our results show 
that the approximation of mono energetic electrons (80 MeV) leads to an over 
estimation of the shower energy by a factor of ~ 8%. 



2 Fluorescence Yield Measurements and Theory 

The fluorescence yield of electrons in air is determined by a competition of two 
factors: the probability of excitation of a nitrogen molecule and the probability 
of its de-excitation via collisions with mainly oxygen molecules (quenching). 
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Based on this assumption Bunner [13] suggested that the fluorescence yield 
(FY) should be proportional to the energy deposit {dE/dx) and can be 
parametrized as: 

dE n 

FYBunner{K, p,T) OC — X 



dx 1 + pBVT 



where K is the kinetic energy of the particle, p is the air density, i? is a 
constant and T is temperature. The dE/dx is the energy loss rate per unit 
path which is dominated by ionization losses. 

More recently, two measurements have been done to determine the fluorescence 
yield as a function of the kinetic energy of the electron. Kamimoto et al. [14] 
measured the fluorescence yield for electrons with energy 1.4, 300, 650 and 
1000 MeV and confirmed the proportionality of the yield with the energy loss 
rate (dE/dx) as suggested by Bunner [13]. Another experiment has been done 
by Nagano et al. [10] and the same proportionality between the yield and the 
dE/dx was verified despite the fact that this experiment was able to make 
measurements at a greater number of wavelength bandwidths than Kakimoto 
et al. did. 

Both measurements suggest that the relation between yield and dE/dx can 
be given by: 



where p is the density in kg m~^ and T is the temperature in Kelvin. The 
constants Ai, A2, Bi, B2 and {dE/dx)K^ is the energy loss rate calculated at 
Kc as given in Table 1. At a given pressure and temperature Kakimoto et al. 
and Nagano et al. differ in the relation between the energy loss rate and the 
fiuorescence yield only by a constant as can be seen in figure 1 where equation 
2 is plotted for T — 300 K and p — 760 mm Hg. 

Taking into account the variation of the density and temperature of the at- 
mosphere as a function of height, it is possible to determine at a given energy 
how the yield varies with height. Figure 2 shows the fluorescence yield as a 
function of height for electrons with kinetic energy of 80 MeV. 

This energy is normally chosen because it is the energy at which the prob- 
ability of an electron undergoing pair-production is equal to the probability 
of absorption by ionization and therefore the number of particles in a shower 
reaches its maximum. However, the energy of the electrons in a shower has a 
wide energy distribution. Furthermore, the energy distribution is different at 
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different fieigfits and for distinct primary particles. 

Nevertheless, the fluorescence yield is a highly dependent function on the 
kinetic energy (through dE/dx) and, since the energy distribution of the elec- 
trons varies with height, the fluorescence yield must vary with height in a 
manner different from that one shown in flgure 2. 

Figure 2 shows that, for 80 MeV electrons, the fluorescence yield is a slowly 
varying function of the height in the atmosphere. The later has been accepted 
as one of the basic features on which rests the fluorescence technique. However, 
as explained above, the characteristic value of the distribution of yields as a 
function of kinetic energy at a given height is not equal to the yield calculated 
at the characteristic value of the energy distribution at a given height. 

In the next sections we are going to convolute the energy distribution of elec- 
trons for heights varying from the sea level up to 7.5 km and for each height 
we are going to calculate the characteristic energy and fluorescence yield. A 
distribution of fluorescence yields for each height was determined allowing us 
to estimate the confldence level of the distribution of yields at each height. 



2. 1 The influence of different atmospheric profiles in the yield 

In applications of the fluorescence yield to measure cosmic rays showers, a 
precise knowledge of the fluorescence yield and its dependences with the kinetic 
energy of the particle are useless if one does not have a good measurement of 
the atmospheric parameters. 

This effect has been shown by B. Keilhauer et al. [15] and we are going to 
stress its importance here by explicitly calculating the ffuorescence yield. 

B. Keilhauer et al. have measured the atmospheric proflle in the Pierre Auger 
Observatory site during an entire year. They have used meteorological radio 
soundings to determine the temperature and pressure variation as a function 
of height. Their results are published in [16] and we are going to use them in 
the following calculation. 

Knowing the pressure and temperature variation as a function of height it is 

possible to calculate the variation of the fluorescence yield as a function of 
height by using equation 2. Figure 3 shows the value of the fluorescence yield 
for three atmospheric models according to Kakimoto et al.. Seasons effects can 
change the fluorescence yield by as much as 4 % [15]. 
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3 The fluorescence yield according to the electron energy spectrum 

In order to evaluate the influence of the distribution of the electron's kinetic 
energy we have done Monte Carlo simulations with the program CORSIKA 
[17] linked to the QGSJET [18] interaction model. We have defined observa- 
tional levels varying from sea level up to 7.5 km and determined the energy 
spectrum of the electrons for each level. 

The simulations were performed with the thinning factor of 10~^ and the 

energy cuts for hadrons, muons, electron and photons were respectively set to 
50 MeV, 50 MeV, 50 keV and 50 keV. Ten vertical showers were simulated for 
each configuration of energy and primary particle. 

Figure 4 shows an example of the kinetic energy distribution electrons of ten 
showers initiated by proton with energy 10^^ eV at 0.5, 3.5 and 7.5 km above 
sea level. It is clear in this figure that the energy spectrum changes significantly 
with height and that this must affect the dependence of the yield with height. 

In order to define a characteristic energy for which we could get a fluorescence 
yield that represents the energy spectrum of electrons we should solve the 
following equation: 

j S{K) X FY{K, p, T)) dK = FY{K,h, P,T) x j S{K) dK (3) 

where in the left hand side we have the energy spectrum {S{K)) convoluted 
with the fluorescence yield {FY{K, p,T)) and in the right hand side we have 
the yield calculated at a characteristic kinetic energy {Kch) times an average 
energy. Using the proportionality of the fluorescence yield with the dE/dx we 
can solve the equation and flnd the values of the characteristic energy which 
we can use in order to calculate a representative fluorescence yield. 

Figure 5 shows an example of the numeric solution for equation 3 for the 
energy spectrum of electrons in proton showers at 1.5 km a.s.l. (above sea level) 
supposing the correspondence of the yield with the energy deposit suggested 
by Kakimoto et. al. 

Equation 3 has two solutions with the same representative yield. We are going 
to quote always the higher kinetic energy. In case of flgure 5 the intersec- 
tion occurs at Kch — 25.8 MeV what leads to a representative yield of 4.36 
photons /m/ electron. 

The same calculation was repeated for all heights, different primary particles 
and energies. 
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3. 1 Varying primary particles type and energy 

For different primary particle types and energies tfie electron energy spectrum 
at a given height may be different and therefore the characteristic energy could 
also change. In order to check the later, we have simulated sets of ten showers 
initiated by proton and iron nuclei with primary energies of 10^^, 10^^ and 
1020 eV. 

Figure 6 shows the electron energy spectrum for proton showers at different 
primary energies at 1.5 km a.s.l.. This figure shows that the number of electron 
at a given energy interval increases with energy independently of the kinetic 
energy. In other words, the overall shape of the spectrum is the same for 
different primary energies what directly implies that the fluorescence yield is 
not a strong function of the energy for primary protons. Exactly the same 
feature happens for primary iron showers where the fluorescence yield does 
not change very much with energy. 

Despite that, small changes in the characteristic energy are seen for different 
primary energies. For example, at 4.5 km for primary proton the character- 
istic energy is 28.6, 31.7 and 32.5 MeV for primary energies of 10^^, 10^^ 
and 10^0 eV, respectively. However, since the dE/dx curve vary slowly with 
energy at this energy range, the correspondent yield changes are within 0.1 
photons/m/electron. The dE/dx values are 4.52, 4.56 and 4.58 respectively. 
Figure 7 shows how the fluorescence yield vary with height for proton primary 
showers with different energies taking into account the energy spectrum as 
explained above for the Kakimoto et al. correspondence between yield and 
dE/dx. We remember that the Nagano et al. correspondence between yield 
and dE/dx is different only a constant factor. 

On average we notice a tendency of the fluorescence yield to get higher with 
increasing primary energy but we emphasize that the differences are always 
smaller than 0.1 photons/m/electron. 

For iron showers wc get the same results and conclusion except that the differ- 
ences between the fluorescence yield for the primary energies are still smaller 
than the ones for proton showers. Figure 8 shows the small differences for the 
yield calculated for primaries protons and iron nuclei with energy 10^^ eV. 

3.2 The fluorescence yield for inclined showers 

The development of a shower is proportional to the amount of matter traversed 
in the atmosphere. Therefore, the energy spectrum of the electrons in a vertical 
shower at 1.5 km is not equal to the electrons energy spectrum in an inclined 
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shower at 1.5 km. The average electrons energy spectrum in showers with 
different inchnations is only equal at the same slant depth. 

Figure 9 shows the fluorescence yield as a function of height for proton showers 
with energy 10^^ eV for three inclinations. 

The calculation regarding the energy spectrum has followed the same proce- 
dure given in section 3. To reproduce the inclination of the shower we took 
the slant depth in the vertical shower and transformed it to equivalent height 
in the inclined shower. This is equivalent to suppose the shower has the same 
energy spectrum at a given slant depth, what would lead to the same charac- 
teristic energy. 

The characteristic energy, as can be seen in equation 3, is independent of pres- 
sure and temperature and can be given as a function of slant depth. However, 
the fluorescence yield calculated for a given characteristic energy is a function 
of pressure and temperature and hence must be calculated at a given height 
for each shower inclination. 

Figure 9 shows that for the same height the fluorescence yield changes as much 
as 0.1 photons/m/electron both from to 37 and from 37 to 53 degrees. 



3.3 Evaluating the fluorescence yield dispersion 



Another important feature of the fluorescence yield of electrons in an air 
shower is their distribution. Historically, experiments have used an average 
of the yield of all electrons at a given height without considering the spread 
of the yield around this average value. 

At each height, we can drawn single electron energies following the energy 
distribution probability as shown in flgure 4. We have drawn 10^ electron 
energies and for each one we calculated the correspondent yield. Figure 10 
shows an example of the distribution of the yield for a proton shower with 
energy 10^^ eV following Kakimoto et al. suggestion. 

As expected, the distribution has a pronounced peak but also has a long tail 
with high yield values calculated for low energy electrons which is normally 
neglected. There is an artiflcial maximum value for the fluorescence yield which 
corresponds to the minimum electron energy in the simulation (50 keV). 
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4 Final Results and Conclusion 



In the reconstruction of cosmic ray showers with the fluorescence technique 
it is impossible to know the energy and the identity of the primary particle 
before the application of the fluorescence yield. Therefore, the safest procedure 
would be to use the average fluorescence yield corresponding to the energies 
and primary particles we expect to measure. 

Figures 11 and 12 show the average fluorescence yield as a function of height 
taking into account the electron energy spectrum. We have averaged over the 
60 showers initiated by proton and iron with primary energies of 10^^, 10^^ and 
10^° eV. The hatched region shows the 68% confldence level of the distribution 
of the yields for each height. 

We also show in the same plot the fluorescence yield calculated for a mono 
energetic electron (80 MeV) for three different atmospheres: US Standard, 
Argentine Winter and Argentine Summer as measured by the Auger Collab- 
oration and published at [16]. 

Only the US Standard atmosphere was used for these calculations. The dif- 
ference from the flxed energy to the spectrum calculations varies with height 
as shown in figure 13 and the average difference is about 8%. 

The exact value the energy reconstructed using the mono energetic model (80 
MeV) has been overestimated varies from shower to shower but our results 
suggests that fluorescence experiments should reevaluate their energy estima- 
tive reducing them by ~ 8%. 

The showcr-to-showcr and particlc-to-particle fluctuation of the yield due to 
the proper consideration of the kinetic energy distribution function of the 
electrons is larger than 10% which at least the double of the variation due to 
atmospheric changes. The distribution of the yield due to the energy distri- 
bution is so big that it would make the measurements of Kakimoto et al. and 
Nagano et al. agree at the 68% C.L. as can be seen in flgure 14. 

The fluorescence yield to be applied to a cosmic ray shower depends on the 
inclination of the axis of the shower. We deflne in equation 3 the concept of 
characteristic energy which is the kinetic energy for which the yield can be 
calculate to take into account the electron energy spectrum in the shower. The 
average characteristic energy for showers initiated by protons and irons with 
energy 10^^, 10^^ and 10^° eV is shown in table 2. The characteristic energy is 
the natural value to be used in the reconstruction of the shower as a function 
of slant depth as it is also shown in flgure 15. 

Given the geometry of the axis of the shower, at a given slant depth the 
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characteristic energy should be used to calculated the fluorescence yield as a 
function of height. Figure 15 shows the points calculated for each height and 
a fit to these points. According to this simple fit the characteristic energy can 
be given as a function of slant depth (%) by the equation: 

Kchix) = 52.028 - 0.04617 x + 1-6968 x 10"^ (4) 

We therefore conclude that the electron kinetic energy spectrum can not be 
neglected in order to calculate the proper fiuorescence yield to be used for 
energy reconstruction of cosmic ray showers otherwise an over estimation of 
8% is artificialy introduced. 
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Constants 


Kakinioto ct al. 


Nagano et ah 


Ai (ni^ kg^i) 


89.0 ± 1.7 


147.4 ± 4.3 


A2 (m2 kg-i) 


55.0 it 2.2 


69.8 ± 12.2 


Bi (m^ kg-i K-°-5) 


1.85 ± 0.04 


2.40 ± 0.18 




6.50 ± 0.3.3 


20.1 ± G.9 


Kc (MeV) 


1.4 


0.85 



Table 1 

Constants used by Kakimoto et al. and Nagano et al. in equation 2. 



Slant Depth (g/cm^) 


Kch (MeV) 


1036 


22.2 ± 0.3 


976 


23.4 ± 0.5 


86 1 


25.1 ± 0.5 


764 


26.7 ± 0.6 


673 


28.3 ± 0.6 


590 


30.6 ± 0.9 


517 


32.1 ± 0.6 


451 


34.7 ± 0.7 


392 


37.5 ± 0.9 



Table 2 

Characterist energy as a function of slant depth. The characteristic energy is the 
average over 60 showers initiated by proton and iron nuclei with energies 10^^ 10^9 
and 10^0 eV. 
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Fig. 1. The dE/dx of an electron in air as a function of its kinetic energy and the 
corresponding fluorescence yield as suggested by Kakimoto et al. and Nagano et al. 
We have used T = 300 K and p = 760 mm Hg. 
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Fig. 2. The fluorescence yield as a function of height in the atmosphere for electrons 
with kinetic energy equal to 80 MeV. Pressure and temperature varies with height 
according to the US Standard Atmospheric Model [16]. 
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Fig. 3. The fluorescence yield as a function of height in the atmosphere for electrons 
with kinetic energy equal to 80 MeV according to Kakimoto et al. suggestion. Pres- 
sure and temperature varies with height according to the models published by B. 
Keilhauer et al. [16]. Arg. Winter and Summer are the average measurements for 
the Pierre Auger Southern site. 
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Fig. 4. Examples of the eleetron kinetic energy distribution for ten proton vertical 
showers with energy 10^^ eV at 0.5, 3.5 and 7.5 km above sea level. 
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Fig. 5. Numeric solution of equation 3 for the electron spectrum for showers initiated 
by protons with 10^^ eV at 1.5 km a.s.l.. 
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Fig. 6. Kinetic energy spectrum for electrons at 1.5 km in ten vertical showers 
initiated by proton with energies 10^^, 10^^ and 10^° eV. 
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Fig. 7. Fluorescence yield according to Kakimoto et al. suggestion as a function of 
height calculated taking into account the energy spectrum as explained above for 
vertical showers initiated by proton with energies 10^^, 10^^ and lO^'^ eV. 
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Fig. 8. Fluorescence yield as a function of height calculated taking into account the 
energy spectrum as explained in section 3 according to Kakimoto et al. suggestion 
for vertical showers initiated by proton and iron with energy 10^^ eV. 
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Fig. 9. Fluorescence yield as a function of height calculated taking into account the 
energy spectrum as explained in section 3.2 according to Kakimoto et al. suggestion 
for showers initiated by proton with energy 10^^ eV and zenith angle of 0, 37 and 
53 degrees. 
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Fig. 10. Distribution of the fluorescence yield calculated taking into account the 
energy spectrum as explained above according to Kakimoto et al. suggestion for 
vertical showers initiated by proton with energy 10^^ eV. 
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Fig. 11. Average fluorescence yield as a function of height. The average was per- 
formed over 30 proton and 30 iron showers with energies 10^^, 10^^ and 10^*^ eV. It 
is shown the yield calculated according to the Kakimoto et al. equation for the US 
Standard, Argentine Winter and Argentine Summer for 80 MeV electron and for 
comparison it is also shown the same calculation for the US Standard taking into 
account the electron energy spectrum. Hatched region corresponds to 68% C.L.. 
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Fig. 12. Average fluorescence yield as a function of height. The average was per- 
formed over 30 proton and 30 iron showers with energies 10^^, 10^^ and 10^° eV. 
It is shown the yield calculated according to the Nagano et al. equation for the US 
Standard, Argentine Winter and Argentine Summer for 80 McV electron and for 
comparison it is also shown the same calculation for the US Standard taking into 
account the electron energy spectrum. Hatched region corresponds to 68% C.L.. 




height (km) 



Fig. 13. Difference between the fluorescence yield calculated according to the energy 
spectrum and calculated for a flxed energy (80 MeV) as a function of height. 
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Fig. 14. Average fluorescence yield as a function of height. It is shown the yield 
calculated according to the Nagano et al. and Kakimoto et al. measurements. We 
have used the US Standard atmosphere and the calculation was done taking into 
account the electron energy spectrum. 




Fig. 15. Average characteristic kinetic energy as a function of slant depth. The 
characteristic energy is the average over 60 showers initiated by proton and iron 
nuclei with energies 10^^, 10^^ and lO^'^ eV. 
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